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Liquid-gas coexistence vs. energy minimization with respect to the density profile in 

the inhomogeneous inner crust of neutron stars 
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We compare two approaches to describe the inner crust of neutron stars: on the one hand, 
the simple coexistence of a liquid (clusters) and a gas phase, and on the other hand, the energy 
minimization with respect to the density proHle, including Coulomb and surface effects. We find 
that the phase-coexistence model gives a reasonable description of the densities in the clusters 
and in the gas, but the precision is not high enough to obtain the correct proton fraction at low 
baryon densities. We also discuss the surface tension and neutron skin obtained within the energy 
minimization. 
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I. INTRODUCTION 

The usual picture of the neutron-star crust [ij is that 
one starts from a Coulomb crystal of nuclei in the outer 
crust. As one goes deeper into the star, the nuclei be¬ 
come more and more neutron rich (as a consequence of 
the increasing electron chemical potential and /? equilib¬ 
rium) until the neutron drip line is reached. This de¬ 
fines the transition to the inner crust, where the nuclei 
(“clusters”) are embedded in a dilute gas of unbound 
neutrons. Descending further into the star, one expects 
to find the so-called “pasta phases” [2l-[l3l|. before one 
finally reaches the neutron-star core where matter be¬ 
comes homogeneous. The inhomogeneous phases in the 
inner crust can also be interpreted in another way Q, 
which is quite common in the study of supernova matter 
(i.e., matter at hnite temperature and out of /3 equilib¬ 
rium) [3 [11 , namely as a consequence of the first-order 
liquid-gas instability of nuclear matter. One aim of the 
present paper is to see to what extent the simple picture 
of phase coexistence can explain certain properties of the 
inhomogeneous phases of the inner crust of neutron stars. 

In particular, the phase coexistence picture corre¬ 
sponds to a hydrostatic equilibrium which can serve as a 
starting point for a hydrodynamic model of the collective 
modes in the inner crust |15l - [l7l| . The collective modes 
can have some impact on thermodynamic and transport 
properties of the crust and may have observable conse¬ 
quences on the cooling of the star [31 ■ While there exist 
several calculations of the collective modes in the crys¬ 
talline phases [l9l - [3| . it seems to be difficult to model 
the collective modes in the pasta phases beyond the hy¬ 
drodynamic approach [I3- 

Within the most simple phase-coexistence picture, 
surface- and Coulomb energies are neglected. An approx¬ 
imate way to include them is the compressible liquid-drop 
model (see, e.g., [ 3 |)- However, this model does not ac¬ 
count for the neutron skin and for the effect of the surface 
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diffuseness on the Coulomb energy. Therefore, we will use 
here a more microscopic approach, namely to parameter¬ 
ize the density profile and determine its parameters by 
minimizing the thermodynamic potential. Nevertheless 
it will turn out that the densities in the dense (“liquid”) 
and dilute (“gas”) regions satisfy quite well the condi¬ 
tions of phase coexistence. 

The article is organized as follows. In Sec. im we de¬ 
scribe the formalisms for the phase coexistence and for 
the more complete minimization of the thermodynamic 
potential. In Sec. mil we show and compare the results 
obtained in both models, and in Sec. nYi we conclude. 


II. FORMALISM 
A. Phase coexistence 


As a hrst approximation, the inner crust could be de¬ 
scribed as a phase coexistence of liquid drops (corre¬ 
sponding to the nuclear clusters) with volume and 
a gas (the dilute neutron gas) with volume ps®-®, which 
satishes mechanical and chemical equilibrium, i.e.. 


pgas _ pliq 

„gas ^ liq 
h^q h^q "i 


( 1 ) 

( 2 ) 


with P® the pressure, and yF the chemical potential of 
neutrons {q = n) or protons [q = p), respectively, in the 
phase i. 

In addition to the neutrons and protons, we consider 
a uniform electron gas to ensure charge neutrality, i.e., 
Vpe = -b where V = -b is the 

total volume and is the proton density in phase i (p®^® 
is not always zero). Instead of working with the volumes, 
it is more convenient to introduce the volume fraction u 
filled by the liquid, which satisfies 

plN Pe - pf^ 

U = —— = —r. - - - . (3) 

^ pI^ - pT 

Furthermore, in neutron stars, matter is in /? equilib- 
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rium, i.e., 


Me = Mn - Mp (4) 

with /ie the chemical potential of the electrons. For prac¬ 
tical purposes, the electron mass can be neglected so that 
Me = (37r^pe)^^^^c. We thus obtain for the volume frac¬ 
tion 


u = 


pfT 


liq 

Pp 


{pn Mp)^ 



(5) 


We use a Skyrme energy-density functional (EDF) (2^ 
to calculate the energy density esk as a function of m„ 
and Pp. This allows us to write explicit expressions for 
the chemical potentials pq = desk/dpq and the pressure 
P = —esk + PnPn + MpMp appearing in Eqs. ([H and ((I|). 
We use different Skyrme parametrizations, all fitted to 
the neutron-matter equation of state: the Saclay-Lyon 
forces SLy4 and SLy7 1^ and the Brussels-Montreal 
forces BSk20 and BSk22 [2^. 


B. Energy minimization 


In the phase-coexistence picture described above, we 
could determine the volume fraction u, but not the ac¬ 
tual size of the clusters. The latter is determined by 
finding the best compromise between the Coulomb en¬ 
ergy (favoring small clusters) and the surface energy (fa¬ 
voring large clusters), which were both neglected in the 
phase-coexistence picture. Let us now consider a more 
sophisticated approach, where these effects are included. 

The surface energy is computed within the semi- 
classical Extended Thomas-Fermi (ETF) approximation 
14 . I^ j. Instead of solving Euler-Lagrange equation 


11|, we use a parametrization of the density profile 
and determine the parameters by minimizing the energy. 
Concerning the Coulomb energy, we use the Wigner-Seitz 
(WS) approximation, i.e., we calculate the Coulomb en¬ 
ergy in an isolated cell r < i?wSi with i?ws chosen such 
that the volume of the WS cell corresponds to the vol¬ 
ume of the unit cell. Inside the WS cell, our density is 
parametrized as 


Pq{r)=pr 


- pT 


l-kexp[(r-rg)/a,] 


( 6 ) 


and the energy is minimized with respect to the 9 pa¬ 
rameters Pq, Pq, Uq, and i?ws- This parametrization looks 
similar to the shape obtained in Hartree-Fock calculation, 
e.g.. Ref. [1^ and has a simple interpretation: p]^^ and 
correspond, respectively, to the asymptotic densities 
in the cluster and in the gas far away from the surface, Vq 
describes the cluster radius, including the possibility of a 
neutron skin if > Pp, and Uq is the surface diffuseness. 
We also tried a more general parametrization using two 
more parameters, namely, raising the Fermi function in 
Eq. dni to a power y,. This additional degree of freedom 


allows one to describe an “asymmetric” surface, but since 
it did not substantially change our results we discarded 
it for simplicity. 

Note that a similar approach was followed by Oya- 
matsu @ and by Pearson et al. [l^]. However, Oyamatsu 
used a completely different parametrization of the den¬ 
sity, which did not become constant inside the cluster. 
Our parametrization (jS]) resembles more the one used by 
Pearson et al. 

It seems very hard to perform the minimization of the 
energy for fixed average densities (averaged over the WS 
cell), since the average densities depend on all the pa¬ 
rameters. We therefore introduce chemical potentials pq 
as Lagrange parameters to fix the densities and minimize 
the thermodynamic potential instead of the energy. To 
be precise, we minimize ut = —P = 0/VwSi where 

fl = Esk ~ PnN — PpZ — PeZ -\- Ef. -|- Ec P E-gx , (7) 

which reduces to 

H = Esk ~ Pn{N P Z) P Ef, P Ec P Egx , (8) 

because of the /3-equilibrium constraint (|T]). In this equa¬ 
tion, Esk denotes the energy obtained with the Skyrme 
functional, N and Z are the total numbers of neutrons 
and protons in the WS cell, Eg denotes the energy of the 
electron gas, Ec the Coulomb energy, and Egx the energy 
due to the Coulomb exchange term. 

The volume Vws of the WS cell depends on Ews and 
on the geometry one considers (spheres, rods, or slabs). 
Let us define the functions Sd{r) and 14(r) describing, 
respectively, the “surface” and the “volume” of a d di¬ 
mensional sphere of radius r, i.e., 

S'i(r)=2, S'2(r) = 27rr, S^{p) = At:p‘^ , (9) 

Vi(r) = 2r, V 2 {p) = Trr^ , V^{p) = ^ttp^ . (10) 

In the case of spheres {d = 3), the volume of the WS 
cell is VVs = f4(Ews)- In the case of rods (d = 2), the 
“volume” Vws = I4(Ews) is actually an area and con¬ 
sequently E, N, Z, etc. represent energies and particle 
numbers per unit length. Similarly, in the case of slabs 
(d = I), the “volume” Vws = ki(i?ws) is a length and 
E, iV, E, etc. represent energies and particle numbers 
per unit area. 

The integrated Skyrme energy in the WS cell, Esk, is 
defined as 

rRws 

Esk=/ d^pn[pn{p),Pp{p)], (II) 
Jo 

with d’^p = Sd{r)dp and H the Skyrme EDF detailed in 
Ref. [ 23 , the kinetic energy density Tq (r) being calculated 

within the ETF approximation as described in [IJ, [ 2 ^ . 
1 


^ Note that terms containing pq pose numerical problems for 
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For the energy of the electron gas, we assume a con¬ 
stant density pe = ^/Vyvs- Hence, the energy can be 
written as Ee = VwsCe- The energy density Ce, taking 
into account the electron mass me, reads [IJ 


Ce 


47r^ 




sinh ^ Xe 


( 12 ) 


with kp^e = (37r^pe)^^^ and Xe = hkp^e/'meC. 

The Coulomb energy Ec is derived from the charge 
density Pc{r) = Pp{r) — Pe- We first calculate the 
Coulomb potential Vc (r) satisfying the Poisson equation 
V^Vc'(^) = iTfe^pc{r) with VVc'(O) = 0. This deter¬ 
mines Vcir) only up to a constant, which is however 
irrelevant for the energy since the total charge of the WS 
cell is zero. The Coulomb energy is given by 


1 rRws 

Ec = -J^ d‘^rpcir)V{r). (13) 

Note that Vcir) and the integral (fOl) must be computed 
numerically, in contrast to the approximate expressions 
one obtains for uniformly charged spheres (or rods, or 
slabs, respectively) with a sharp surface which are often 
used in the literature iiiHl- 

We relaxed the assumption of a constant electron den¬ 
sity by using a screened Coulomb potential instead of the 
full one, i.e., replacing in the Poisson equation Vcir) 
by (V^ —l/A^)Vc'(r), where A denotes the Debye screen¬ 
ing length, 1/A^ = (4Q;/7r)fc|,g (neglecting the electron 
mass), with a = jhc. However, as already noticed in 
, the screening length of the electrons is so large that 
this effect is negligible and the calculations of the present 
paper were done without electron screening. 

The exchange term of the Coulomb interaction is com¬ 
puted within the Slater approximation [TTj | . For the pro¬ 
tons, it is given by the integral of 


eex.p(r) = e^pt^^ir), (14) 

over the WS cell. For the relativistic electrons, the ex¬ 
change energy includes also contributions from transverse 
photons a nd g ets positive for Xe ^ 2.53 [2^. Its expres¬ 
sion reads ISOl 



(15) 


d = 2 and 3 since the parametrization of the density 10 is only 
approximately constant at the origin (in contrast to the density 
profile used in Ref. [^ . whose derivatives vanish at r = 0 and 
r = Rws)‘ We avoid these problems by integrating the terms 
of the functional containing pq by parts and then using the 
form of the functional obtained in this way in the numerical 
calculation. 
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FIG. 1; (Color online) Neutron and proton densities in the 
coexisting liquid (droplets) and gas phases, calculated with 
the Skyrme parametrization BSk22. 



FIG. 2: (Color online) Same as Fig. [T]but with a logarithmic 
scale for the baryon density. 


III. RESULTS 
A. Phase coexistence 

Let us hrst discuss the results obtained within the sim¬ 
ple phase coexistence approach. Figured] represents the 
densities in the gas and in the liquid as functions of the 
baryon density 

p, = uip^^ + pD'^) + (I - u)ipr + pT) > (16) 

obtained with the BSk22 interaction Results cal¬ 
culated with other Skyrme parametrizations are close to 
those displayed for BSk22 (cf. Fig. [5| for SLy4). The 
main difference is the transition to uniform matter, i.e., 
the point where the liquid fills the whole volume (u = I). 
For instance, the transition occurs at pb « 0.09 fm“^ 
with SLy4, while with BSk22 it is below 0.08 fm“^. 

Figured) shows the same results, but on a logarithmic 
scale so that the low density region is better visible. The 
region below pb « 10“^ fm“^ corresponds to the outer 
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FIG. 3; (Color online) Proton density in the gas phase (and 
later in uniform matter) as function of the total baryon den¬ 
sity calculated within the phase coexistence approach using 
the SLy4, SLy7, BSk20 and BSk22 Skyrme parametrizations. 

crust, where = 0 (i.e., p,n,Pp < 0). In the 

limit of extremely low densities, the droplets are made 
of symmetric nuclear matter at saturation density (p„ = 
Pp = 0.08 fm“^). The kink in the densities visible in 
Fig. [5] corresponds to the transition between the outer 
and the inner crust, i.e., to the point where becomes 
positive and the neutron gas appears. 

Close to the transition to uniform matter, also p,p be¬ 
comes positive and the gas phase does not only contain 
neutrons, but also a small amount of protons. Figure [3] 
zooms on densities where protons are present in the gas. 
Note that the proton density in the gas is always very 
low, less than 10“^ fm~^ for all Skyrme parametriza¬ 
tions, compared to the other densities calculated in the 
phase coexistence framework. 

B. Energy minimization 

In the preceding subsection, we treated the inner crust 
as two nuclear fluids in phase coexistence. As already 
mentioned in Sec. IIIBl this approach misses surface and 
Coulomb effects, which are included in the minimization 
of the thermodynamic potential with respect to the pa¬ 
rameters of the density profile given in Eq. (|6]). Actually, 
we perform the minimization for each of the three dif¬ 
ferent geometries discussed in Sec. IIIBl namely spheres 
(3D), rods (2D), and slabs (ID). The geometry giving the 
lowest thermodynamic potential should be the one that 
is physically realized. 

From now on, we restrict ourselves to the range > 0 
corresponding to the inner crust and show the results as 
functions of the average baryon density pb = A/VvvSj 
where A = N + Z \s the total number of nucleons in the 
WS cell. Figure m displays the difference w — between 
the thermodynamic potentials uj in 3D, 2D, and ID geom¬ 
etry, obtained by numerical minimization with the SLy4 
interaction, and the thermodynamic potential Uu of uni¬ 
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FIG. 4: (Color online) The thermodynamic potential lo as 
function of the baryon density pb, with the thermodynamic 
potential of uniform matter subtracted. In the inset, a phe¬ 
nomenological function aJ is subtracted, which approximates 
the average behavior of the three tj’s. The SLy4 interaction 
was used in this calculation. 


form npe matter in /3 equilibrium with the same baryon 
density. The difference is negative, confirming that the 
inhomogeneous phase is favored over uniform matter in 
this density range. In order to make the differences be¬ 
tween the 3D, 2D, and ID geometries better visible, we 
subtract in the inset a purely phenomenological func¬ 
tion UJ which approximates the average behavior of the 
three w’s. At densities below 0.061 fm~^, the most favor¬ 
able phase is the crystal (3D). From 0.061 to 0.073 fm~^ 
the preferred phase are the rods (“spaghetti”, 2D). Fi¬ 
nally between 0.073 and 0.081 fm“^ we find the slabs 
(“lasagne”, ID), until the system transforms into uniform 
matter. In contrast to other work Hi, [mil , we did 
not find “inverted” geometries such as tubes and bubbles 
(“Swiss cheese”), which would correspond to p^^^ > p^^’^ 
in the parametrization (jH]) of the density profile in 2D 
and 3D, respectively. 

Let us note that the energy differences between the 
three geometries are extremely small compared to the 
total energy, especially between 2D and 3D beyond ^ 
0.05 fm“^, so that one may expect coexistence of dif¬ 
ferent geometries. This is because, in contrast to most 
studies of the pasta phases, we do not consider a fixed 
proton fraction, but /? equilibrium. As pointed out in 
Ref. im, the small proton fraction (see below) corre¬ 
sponding to P equilibrium is very unfavorable for the 
formation of pasta phases. Our method results necessar¬ 
ily in first-order phase transitions between the different 
geometries. In reality, however, it might happen that the 
system passes continuously from one phase to another, 
e.g., by deforming the nuclei in the 3D phase before they 
merge into rods H HI . 

The densities p* in the favored geometry obtained by 
the minimization are displayed in Fig.[5]as the solid lines. 
We notice a discontinuity at the transition from 2D to 
ID. It is dne to the first-order phase transition men- 
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FIG. 5: (Color online) Solid lines: densities in the cluster 
and in the gas corresponding to the geometry 

(3D, 2D, ID, or uniform matter) that minimizes cj in Fig. 2] 
obtained by the minimization of the SLy4 functional. Dots: 
results obtained within the phase coexistence approach as in 
Fig. m but with the SLy4 interaction. 

tioned above (the corresponding jump between 3D to 2D 
is too small to be seen on the figure). For comparison, 
the dotted lines in Fig. [5] were calculated as in Fig. [1] 
within the phase coexistence framework. The densities 
obtained within both approaches are quite similar, which 
means that even with Coulomb and surface effects, the 
mechanical and chemical equilibrium, Eqs. o and ©, 
are approximately satisfied for the densities in the clus¬ 
ter and in the gas away from the interface. The main 
difference is the transition density from the inner crust 
to uniform matter, i.e., to the neutron star core. Since 
Coulomb and surface effects favor uniform matter, the 
transition happens earlier (i.e., at lower pb) in the mini¬ 
mization than in the phase coexistence approach. 

Although the differences in the densities are small, they 
have a sizable effect on the proton fraction Yp = ZjA. 
This can be seen in Fig. [6j where we compare the pro¬ 
ton fractions obtained within the phase coexistence ap¬ 
proach (solid lines) and by energy minimization (dots). 
It appears that at low baryon density, the proton fraction 
obtained by minimizing the energy is significantly lower 
than in the phase coexistence picture, although the pro¬ 
ton densities in the liquid are very close (cf. Fig. [S]). A 
very similar difference between the two approaches is ob¬ 
served in the Relativistic Mean-Field (RMF) framework 
0. This disagreement can be traced back to the tiny 
difference in due to the small Coulomb and surface 
corrections to the chemical and mechanical equilibrium. 
Since in this region the total density is dominated by the 
density of the gas, pb ~ p^, and the volume fraction u 
is approximately proportional to the difference ps — 
the quantities u and consequently also Yp are very sensi¬ 
tive to small deviations of pf^^. 

The proton fractions obtained by energy minimization 
are similar to those obtained by Pearson et al. within 
the Hartree-Fock-Bogoliubov (HFB) model with BSk in¬ 


FIG. 6: (Color online) Proton fraction Yp as function of the 
baryon density pb. The proton fraction is computed using the 
SLy4 interaction with the phase coexistence (solid line) and 
the energy minimization (dots) approach. 


teractions. Bowever, in contrast to the BFB model, the 
ETF approximation does not include shell effects. There¬ 
fore, our proton number in the 3D phase varies smoothly 
from Z ss 34 at low baryon density to Z Ri 26 at the 
transition to the 2D phase, while the BFB model gives 
Z = 40 in the whole inner crust except in the case of 
the BSk22 interaction where Z jumps from 40 to 20 at 
Pb = 0.035 fm-3 m. 

C. Properties of the liquid-gas interface 

Let us discuss in some more detail the properties of the 
liquid-gas interface obtained by energy minimization. To 
that end, we have to compare our WS cell with a ref¬ 
erence system containing a cluster with constant densi¬ 
ties and a sharp surface, surrounded by a gas with 
constant densities The presence of a neutron skin 

Sn = fn — Tp > 0 complicates this comparison, and we fol¬ 
low Douchin and Baensel who discussed this problem 
in detail. 

We define the radius r®® of the reference cluster in such 
a way that it contains the same number of protons as the 
actual WS cell, i.e., 

(Pp" - Pr)Vd{rf) + pfVws = ^ (17) 

{Pp^^ = 0 in most cases, except near the crust-core tran¬ 
sition, cf. Fig. ED. Note that r®® coincides with Tp in the 
case d = 1, but not in the cases d = 2 or 3. If we used 
an asymmetric surface (see discussion below Eq. (O), r®® 
would differ from also in the case d = 1. 

Since we define the reference cluster with a common 
surface at r®® for protons and neutrons, i.e., without 
neutron skin, the reference system contains less neutrons 
than the actual WS cell. Therefore, rather than compar¬ 
ing the energies of two systems having different numbers 
of particles, one should compare their thermodynamic 
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FIG. 7: (Color online) The surface tension cr as a function of 
the baryon density ph. For each pb, the result corresponding 
to the most favored geometry is displayed. 
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FIG. 8: (Color online) The number of the additional neutrons 
due to the skin, Ns, divided by the cluster surface Sd(r'^)- 


potentials. The surface contribution to the thermody¬ 
namic potential, VLg, is defined as the change in (ex¬ 
cluding Coulomb) with respect to the “bulk” thermody¬ 
namic potential (i.e., excluding gradient terms) of the 
reference cluster. We denote = 14(r®®) the volume 
of the reference cluster and = Vy^g — I/i't the vol¬ 
ume of the surrounding gas. Then the surface potential 
can be written as 


- - Mnpr - win > (i*) 

where egj^ is the energy density obtained with the Skyrme 
functional Ti in the case of uniform matter with densities 
Pn and pI- 

Analogously to Eq. (fT71) , we define an effective neutron 
radius r®®, and an effective volume of the neutron skin 
Vs = Vd(r®®) — Vd(r®®). Then the number of neutrons in 
the skin. Ns, is given by 

IV« = E.(p^-pr), (19) 

and Eq. m can be rewritten as 

Qs = Esk - - ys-ggas _ ^ 


Finally, the surface energy Es is given by Eg = Qs+PnNg 


The surface tension is approximately given by cr = 
lls/«5'(i(r®®), where Sd{r) is defined in Eq. ([HI). Although 
Vd{r) and Sd{r) have the dimensions of a volume and 
of an area only in the case d = 3, the ratio ns/S'd(r®®) 
is always an energy per area. Note that the above def¬ 
inition of the surface tension is only exact in ID or if 
the cluster radius is so large that curvature effects can 
be neglected. In Fig. [7] we display the surface tension. 
We see that it decreases with increasing pf,. This is not 
surprising because with increasing pb the densities in the 


gas and in the liquid get closer to each other. In the 
low-density limit, where we have essentially isolated nu¬ 
clei, our surface tension is consistent with the surface 
energy in the Bethe-Weizsacker semi-empirical mass for¬ 
mula from which one obtains a = 1.03 MeV fm“^ [s^. 
Our surface tension is similar to the RMF results shown 
in Fig. 6(c) of Q. 

Another interesting quantity is the number of addi¬ 
tional neutrons in the skin. Again, we normalize it to the 
surface and display in Fig.[5]the number of skin neutrons 
per unit surface, Ng/Sd{rp^)- We observe that, after a 
maximum at pb ^ 0.014 fm~^, the number of neutrons 
in the skin decreases with increasing baryon density, al¬ 
though the skin thickness remains about ~ 0.45 fm, 
because the density difference — p®^® decreases. Note 
that, compared to the total number of neutrons in the 
WS cell, the number of neutrons in the skin is very small: 
for instance, for pb = 0.02 fm“^ there are about 800 neu¬ 
trons in the cell, 200 in the cluster, and only 23 in the 
skin. 


IV. CONCLUSIONS 

In our work we compared the results obtained for 
the inner crust within a simple phase-coexistence pic¬ 
ture with results of the minimization of the energy with 
respect to the density profile. We used Skyrme interac¬ 
tions and the ETF approximation to calculate the energy. 
Both approaches give similar results for the neutron and 
proton densities in the clusters and in the gas, despite the 
Coulomb and surface effects included in the minimization 
but not in the phase coexistence. 

However the phase-coexistence picture is insufficient 
to predict the cluster sizes and transitions between dif¬ 
ferent geometries, since these result from the competi¬ 
tion of Coulomb and surface energies not included in the 
phase coexistence. With increasing baryon density, we 
find crystals, rods and plates, but no “inverted” geome- 
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tries such as tubes and bubbles. Because of the small pro¬ 
ton fraction in /3 equilibrium, the energy differences be¬ 
tween the different geometries are extremely small. An¬ 
other effect of Coulomb and surface energies is to shift 
the crust-core transition to lower baryon density. 

Although the densities obtained by energy minimiza¬ 
tion do not present strong deviations from those of the 
phase coexistence, the proton fraction does, especially at 
low baryon densities. The small corrections to mechan¬ 
ical and chemical equilibrium due to Coulomb and sur¬ 
face effects slightly modify the density of the neutron gas, 


resulting in a considerable reduction of volume and pro¬ 
ton fractions. The proton fractions obtained within the 
minimization are similar to HFB results in the literature 
[s^ . although shell effects, which are present in HFB, are 
missing in the ETF approximation to the energy. 

The minimization allowed us to calculate the surface 
tension and the number of skin neutrons in the cluster 
surface. We plan to use these results in an improved 
hydrodynamic description of collective modes in the in¬ 
ner crust similar to the one of Ref. da but including 
Coulomb and surface effects. 
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